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ABSTRACT 

Flaw detection and characterization with thermographic techniques in graphite polymer composites is often 
limited by localized variations in the thermographic response. Variations in properties such as acceptable porosity, 
variations in fiber volume content and surface polymer thickness result in variations in the thermal response that 
in general cause significant variations in the initial thermal response. These variations result in a noise floor 
that increases the difficulty of detecting and characterizing deeper flaws. The paper investigates comparing 
thermographic responses taken before and after a change in state in a composite to improve the detection 
of subsurface flaws. A method is presented for registration of the responses before finding the difference. A 
significant improvement in the detectability is achieved by comparing the differences in response. Examples of 
changes in state due to application of a load and impact are presented. 
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1. INTRODUCTION 

With the increased application of composites in commercial and military aircraft, rapid large area inspections are 
more important for ensuring the safety and reliability of aircraft. Graphite fiber reinforced composite materials 
are being used as primary structure due to their high stiffness and strength to weight ratio. One particular 
interest is the detection of delaminations that can appreciably reduce the compressive strength of a composite. 

Thermography is a demonstrated technique for rapid inspections of materials and structures and has been shown 
to have great potential for detection of delaminations in composites . 1-6 These efforts have included a variety 
of heating and data reduction techniques to improve the detectability and assessment of the size and depth of 
delaminations. 

Single sided thermographic detection of deeper delaminations in graphite fiber reinforced polymer(GFRP) com- 
posites is difficult . 7 There are three obvious reasons for this. The first reason is GFRP composites are an 
anisotropic material, with the fibers typically in the plane parallel to the surface. The high thermal conductivity 
of the graphite fibers relative to the thermal conductivity of the polymer results in the heat diffusing much faster 
in the direction of the fibers. The in-plane diffusion is much faster than the through the thickness diffusion, 
therefore the heat diffuses around deep flaws and does not result in a significant temperature increase at the 
surface. A second reason deep delaminations are difficult to detect is the inherent inhomogenetity of the GFRP 
composite. Small pockets of trapped air near the surface or variations in the fiber volume content results in 
variations in the surface temperature that can be as large as the variation caused by a delamination. Thirdly for 
depths close to the thickness of the composite, the maximum possible difference between the delaminated and 
undelaminated composite is proportional to the ratio of the thickness of the composite below the delamination 
to the total thickness above the delamination. As this thickness below the delamination goes to zero, so does 
the difference. 

A fourth, perhaps less obvious reason deep delaminations are difficult to detect in composites is a result of the 
nature of delamination in the composites. Delaminations in composites are a thin air gap between two plys in 
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the composite. If the thickness of the air gap is small, the heat flow across the gap may not provide a significant 
enough thermal resistance relative to the layer above the delamination. The effect of the thermal resistance on 
the thermal response of the delamination has been previously examine by several groups. 8-10 In Sec. 2, it is 
shown that the contact resistant of the delamination must be significant relative to the thermal resistance of the 
layer above it to be detectable. 

This paper investigates the potential for increasing the detectability of flaws in composites by performing a 
“difference thermography.” For difference thermography, the thermal response of the composite is measured 
while the composite is in one state, referred to as a “baseline” measurement. The state of the composite is then 
changed and by carefully aligning the images corresponding to the thermal response with an affine transform, a 
point-by-point comparison of the thermal responses is made. By performing the affine transform on the data, 
the specimen and measurement system are not required to stay in fixed positions relative to each other for the 
baseline and post change-in- state measurements as was done previously. 11 Two examples are given, the first 
using a wedge to open an existing delamination to change the state of the composite. In the second example, 
the change in state is pre and post impact. For both cases, calculating the point-by-point state change induced 
difference in thermal response is shown to significantly improve the detectability of subsurface damage. 


2. ONE DIMENSIONAL MODEL FOR COMPOSITE WITH DELAMINATION 


A simple analytic solution does not exist for the one-dimensional heat flow in a multilayered material. The Laplace 
transform solution for two layers of thickness L\ and L 2 coupled by an intermediate contact resistance (R) is 
known. Since the configuration of interest is a composite with a delamination, the first and second layers are 
assumed to have the same thermal conductivity (. K ) and diffusivity (a). For the surface with the incident heating, 
the Laplace transform of the temperature response of the front surface (Tp(s)) for instantaneous heating is given 
by 12 
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Tb(s ) is the Laplace transform of the back surface temperature and / is the incident flux. Solving for Tp(s) 
gives 
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The Laplace transform for the difference between delaminated and undelaminated thermal responses(i? = 0) is 
given by 
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Note that the full thickness response is the reference being subtracted and not the the Laplace transform of a 
semi-infinitely thick thermal response media to instantaneous heating or f^/a/s/K. 

A simple analytical solution for the inverse Laplace transform of Eq. 3 does not exist, however if one assumes 
the contact resistance is small, then the time response is given by 
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Figure 1, Time dependence D(t ) as expressed in Eq. 9 which gives relative change in the thermal response 
of a composite due to a delamination if the contact resistance of the delamination is small relative to the 
thermal resistance of the layer above it. For comparison, the full series solution (G(t,p)) for different values 
of p are also plotted. 
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and L\ — pL and L2 = L(1 — p) with L being the full thickness of the composite. The first two terms inside 
the bracket correspond to the thermal response of a layer with no delamination. The third term is the difference 
between the single layer response and the response with a delamination. This term must be significant, relative 
to 1 for the delamination to produce a significant change in the thermal response of the composite. If p < 0.3 
(delamination is in top half of the specimen) , to a good approximation, this third term can be further simplified 
to 
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where r = a/ L\. The time dependence of D(t) is shown in Fig.l. D(t) has a maximum value at t = r of e -1 times 
KR/L i, which is the ratio of thermal resistance of the layer above the delamination to the contact resistance. 
After reaching this maximum, since the contact resistance is small, the heat flows through the delamination. 
This reduces the front surface temperature to the value for no delamination without requiring any heat flow 
around the delamination. For deeper delaminations, the effect of the delamination is reduced since L\ is greater 
and G(t,p) has a lower maximum value. 


3. METHOD FOR REGISTRATION OF TWO THERMAL RESPONSES 

The two thermal responses can be considered to be a series of images, Ai(x, y ) and Bi(x, y ), where i corresponds 
to time of the image and x and y correspond to the location of pixels in the image. To register the two thermal 
data sets, the first unsaturated thermal images (defined as i = 1 for each thermal response are registered to 
each other. For the cases examined, the baseline data sets are held fixed and the data sets from post change are 
transformed for registration with the initial state. For many cases a simple rotation and translation is required. 
However, it is possible that between the data acquisition, the configuration could have changed enough that the 
plane of the specimen is in a plane that is rotated relative to initial configuration. For those cases, an affine 
transform is required. 



Figure 2. Registration of images based on affine transform, (a) Infrared image of specimen tilted out of the 
normal plane of data acquisition by approximately 30° and rotated by approximately 16°. (b) Results of the 
affine transform using T. (c) Fixed reference infrared image with tilt ~ 0 and rotation « 0 that was the 
target of the the optimization routine. 


The affine transform is given by 
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where x and y are the coordinates of the initial frame of reference and x' and y' are the coordinates of the 
transformed frame of reference. The elements of the matrix, (ai ? i • • • 613,3) are seven independent parameters 
which represent the affine transform. 


When the transformation is a simple rotation by 0 , followed by a translation of x t and y t and finally a scaling 
then Eq .10 becomes 

cos( 0 ) sin{ 0 ) x t 
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where m is the magnification. 

Registration of the two thermal responses is performed by selecting a region on interests in Ai(x , y). To determine 
the proper value for ip = [ai,i, 621,2, 621,3, a 2,i, <^2,2 , 6*2,3, 613,3], initial values are chosen for the different elements 
of ip and image Bi(x,y) is transformed to Ci(x,y). The pixels of Ci(x,y) that correspond to the pixels of the 
region of interest in Ai(x,y) are amplitude and offset matched using a least squares estimation. The sum of 
the squared differences of the least squared estimation is used as the cost for a simulated annealing routine that 
varies ip to determine T, the value of the vector corresponding to the global minimum for the cost. The T is 
then used as the parameters for performing the affine transform of Bi(x,y) for i = 1 to N (number of images) 
to Ci(x,y). Ai(x,y) is then subtracted from Ci(x,y) for i = 1 to N to calculated the difference thermography 
data set. 
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An example of the results of this process is shown in Fig. 2 . A composite specimen with a wedge insert into a 
delamination was tilted out of the typical measurement plane by approximately 30 ° and rotated by approximately 
16 °. This should be considered to be an undesirable initial alignment, however, it is presented as a demonstration 
of the capability of the registration technique. As can be seen in Fig. 2 , the registration technique results in an 
image that has excellent image registration with the target image. The technique resulted in excellent registration 
for all cases considered. 


4. APPLICATION OF DIFFERENCE THERMOGRAPHY TO COMPOSITE 

SPECIMENS 

Difference thermography was applied to two different types of composite damage. The first type of damage 
examined was an existing delamination with and without a wedge insert into the delamination to increase the 
gap width. The second type of damage was impact damage, where the difference was taken between pre and 
post impact data. 




Estimated Gap Width-0.1 8 mm 



Figure 3. Side view of the composite specimen discussed in Sec. 4.1. (a) The composite specimen with 
existing delamination, (b) Same specimen with wedge inserted into the delamination. The estimated width 
of the air gap is shown at the same point of the specimen. 




(b) 

Figure 4. Difference between thermal response of composite with wedge opening delamination and reference, 
(a) Reference is the mean of the thermal response of in a region without delamination, (b) Reference is 
point in the aligned thermal thermal response, taken before delamination is widen with wedge. For both (a) 
and (b), the difference signals acquired from the undamaged region of the specimen is shown in blue and 
difference signals from the delaminated region is shown in red. 





4.1 Composite with delamination 

The first case examined was a composite specimen with an existing delamination. A picture of the side of the 
composite with and without the wedge inserted is shown in Fig. 3. As can be seen from the figure, the width of 
the air gap decreases as one moves from the left side of the image toward the middle of the specimen. At the 
points with the estimated gap widths, RK/Li « 3 for Fig. 3(a) and RK/Li & 8 for Fig. 3(b). 

After performing the image registration as discussed in Sec. 3, the point-by-point differences in the thermal 
response were calculated. A comparison of two methods for calculating the difference in thermal response are 
shown in Fig. 4. For both Fig. 4(a) and Fig. 4(b), the data shown are from the composite with a wedge inserted 
into the delamination. Fig. 4(a) displays the difference between thermal responses and a reference calculated 
from the mean value at each time in a region of the composite where no delamination exists. The red response 
difference signals are from the delaminated region of the composite and the blue difference signals are from the 
undamaged region of the composite. As can be seen from the figure, even with the significant increase in gap 
width, there is considerable overlap between the difference signals of delaminated and undelaminated regions for 
all times. 

Instead of using a single reference single, the difference signals were calculated from the signals of the registered 
data sets from prior to and post insertion of the wedge into the delamination. The results are shown in Fig. 4(b). 
By differencing baseline and registered post thermal responses, the variability of the undelaminated region of the 
specimen is significantly reduced. This results in a difference response for the delaminated region that is greater 
than the variability of the undelaminated region from approximately 5 seconds to the end of the time record at 
10 seconds. 

Averages of the thermal responses of the composite from 3 to 5 seconds with and without the wedge in the 
delamination are shown in Fig. 5 (a) and Fig. 5(b) respectively. When there is no wedge inserted into the delami- 
nation, it is difficult to detect the delamination in the averaged thermal response. After inserting the wedge, the 
delamination is more obvious in the image, however, it is difficult to establish where the edge of the delamination 
is. There is a gradual change in the response as one moves from the undelaminated to delaminated region of the 
composite. 

The averaged response obtained by calculating the difference of the registered thermal responses for prior to 
and post insertion of the delamination is shown in Fig. 5(c). As can be seen from this figure, the differencing of 
the thermal responses considerable improves the visualization of the delaminated region. It is also much more 
evident where the delamination ends on the specimen. 

If one assumes the baseline thermal response of the composite is approximately the thermal response of the 
undelaminated composite, then one can use Eq. 6 to estimate the depth of the delamination from the difference 
response. The inversion of Eq. 6 is performed numerically using the Talbot method. 13 One of the unknowns in 
Eq. 6 is i?, however for large gap t hicknesses (~0. 01 cm or RK/Li ~ 2 , the shape of the thermal response is not 
very dependent on the value of R for times less than r ^ 20 second, therefore RK/Li was set to 2. The good 
agreement of model and data is shown in Fig. 6. From fits of a significant portion of the delaminated region, the 
depth of delamination from the thermal fit was estimated to be 0.22 d= 0.02 cm which is in agreement with the 
depth of delamination from the measurement of the specimen of 0.219 cm. 

4.2 Composite with impact damage 

A second case examined was a composite structure undergoing impact testing. Baseline and post impact data 
were acquired on the structure. The first unsaturated images of the two data sets were registered using the 
method discussed in Sec. 3, then the point-by-point differences in the thermal response were calculated. There 
were reference marks on the surface of the composite, however they were not required to register the data sets, 
since the initial images had significant variations due to local anomalies in the composite. The markings were 
intended to indicate where the composite was to be impacted. 

For one of the impacted regions, the thermal indications of impact damage for this case are not significantly 
larger than the normal variations in the composite as can be seen by comparing Fig. 7(a) and Fig. 7(b). When 
the difference thermography image is examined(Fig. 7(c)), the impact damage becomes the predominate feature 
in the image. What also becomes obvious is damage that is not visible without performing the difference 




Figure 5. Average from 3 to 5 seconds of singled sided thermal response, (a) Results of the averaged signal 
for composite without wedge in delamination, (b) Results of the averaged signal for composite with wedge 
in delamination, (c) Results of the difference between the thermal response of the same specimen with and 
without the delamination gap increased. Arrow indicates the approximate end of the delamination, based 
on examination of the edge of the composite. 
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Figure 6. The difference between thermal responses with and without the delamination gap widened for 
point over delamination and point over no delamination. The fit of the delamination data based on Eq. 6 is 
overlaid on the delamination response. 



(a) (b) (c) 

Figure 7. Results of difference thermography technique for impacted specimen, (a) Infrared image of region 
of structure taken prior to impact, (b) Infrared image of same region taken after impact, (c) Difference 
image of region at 0.5 second following flash heating taken prior to impact. All images represent data at 0.5 
second following flash heating. 


thermography. The depth of the damage was not determined from this data, due to the small lateral size of the 
damage. 


5. SUMMARY 

Difference thermography is shown to be an effective method for improving the detectability of flaws in a composite. 
Two examples are discussed, the first using a wedge to open an existing delamination to change the state of the 
composite and the second the change is state is pre and post impact. For both cases, the difference thermography 
is shown to significantly improve the detectability of subsurface damage. An affine transform of infrared images 
enabled registration of the thermal response without requiring the specimen and measurement system stay in 
fixed positions relative to each other for the prior to and post change in state measurements. 
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